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ABSTRACT 

We present a time-dependent wave-packet propogation calculation for the H 2 + CH Z -> 
H + CH 4 reaction in six degrees of freedom and for zero total angular momentum. Initial 
state selected reaction probability for different initial rotational-vibrational states are pre- 
sented in this study. The cumulative reaction probability (CRP) is obtained by summing 
over initial-state-selected reaction probability. The energy-shift approximation to account 
for the contribution of degrees of freedom missing in the 6D calculation is employed to ob- 
tain an approximate full-dimensional CRP. Thermal rate constant is compared with different 
experiment results. 
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I. INTRODUCTION 




The reaction H + CH A —¥ Ho + CH 3 and its reverse reaction CH 3 + H 2 -4 H + CH 4 
are of obvious importance in methane pyrolysis and combustion [1]. In addtion, the reverse 
reaction also plays an important role in the modeling of the chemistry of hydrocarbons 
in the atmospheres of giant planets [2], Numerous experimental [3-12] and theoretical 
studies [13-18] have been reported on these two reactions. Recently, the forward reaction 
H + CH 4 — > H 2 + CH 3 has become the focus of a number of quantum dynamics calculations 
[19-26] in reduced dimensionality and a full-dimensional transit state flux method quantum 
calculation. Nevertheless, no quantum dynamics studies have been reported on the reverse 
reaction H 2 + CH ?J — » H + CH 4 so far. 

In the reduced dimensional [19,20,23,24] calculations for the forward reaction H+CH 4 — » 
H 2 + CH 3 , usually energy shifting [27-29] methods are needed to obtain approximate full di- 
mentionality cumulative reaction probability for zero angular momentum J, and J-K shifting 
methods [27-29] are employed to calculate total cumulative reation probability for the full 
dimensionality rate constants for both reduced dimensional and full dimensional calculations 
[22,25]. In the reduced dimensional quantum calculations, initial stated selected reaction 
probabilities were calculated. However, for the full dimensional calculation, the calculations 
were limited to the CRP only due to the limition of the transit state flux method calcula- 
tion to obtain initial state selected reaction probability, as well as due to the difficulty of 
calculations in full dimensionality. 

All the quantum calculations mentioned above were done using the twelve degrees of 
freedom potentional energy surface of Jordan and Gilbert [30]. This potential has a C 3v 
transition state with a collinear structure of H — H - C and the center of mass of H 3 . The 
vibrationally adiabatic barrier height of this potential energy surface is 10.9 kcal/mol for 
the forward reaction, H -I- CH 4 — * H 2 + CH 3 and 11.0 kcal/mod for the reverse reaction, 
H 2 + CH 3 — > H + CH 4 , respectively. This potential energy surface is based in part on 
previous ab initio calculations of the saddle point, as well as on a previous, semiglobal 
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potential used in variational transition state theory calculations fo the rate constant. The 
rate constant from the theoretical calculations of the forward reaction, comparing to the 
experiment, indicates that the barrier height on the Jordan-Gilbert potential may be too 
low. 

In this paper, we report six degrees of freedom, time -dependent quantum calculations 
of the reaction H 2 T CH 2 — > H -I- CH 4 , by treating the reaction as a diatom-diatom reac- 
tion system and carry out a full dimensional quantum dynamics calculation for this pseudo 
diatom-diatom reaction. Initial state selected reaction probability, cumulative reaction prob- 
ability(CRP), and thermal rate constant are calculated. We apply energy shifting method 
to obtain the thermal rate constants and compare our results with experiments. 

This paper is organized as follows: Section II briefly reviews the theory used in this 
calculation and the energy-shifting methods are also presented. Section III gives numerical 
results of the calculation, and section IV concludes. 


II. THEORY AND METHODS 


A. Time-dependent wave-packet propogation method 


This study employs the time-dependent wave-packet propogation method [31-33] to carry 
out the quantum dynamics calculation. The reverse reaction, H 2 + CH 3 —¥ H + CH 4 , is 
treated as a diatom-diatom reaction, H 2 + CX — > H + HCX, where X is a pseudoatom with 
the mass of the three unreactive hydrogen atoms and located at their center of mass. The 
vibrational motion of CX approximately represents the umbrella motion of CH 2 , where the 
three hydrogen atoms move in phase , and the C — H bond length is fixed at its equalibrium 
geometries. Thus the reduced dimensionality calculation for this pseudo diatom-diatom is 
six degrees of freedom. 

The Hamiltonian for this model system in the reactant Jacobi coordinates is given by 


ti 2 d 2 ( J - j ) 5 
2 fj,dR 2 + 2 nR 2 


+ ^ri ( r i) + K 2 {r 2 ) + 


V(R,r 1 ,r 2 ,6i,9 2 ,<f>), 
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where p, is the reduced mass of H 2 with respect to the center of mass of the pseudo diatom 
CX , ri is the distance for H 2 , r 2 is the distance between C and X, R is the distance from the 
center of mass of H 2 to the center of mass of CX, J is the total angular momentum operator, 
jx and j 2 are the rotational angular momentum operators for H 2 and pseudo diatom CX 
respectively, 9\ and 0 2 are the two Jacobi angles between R and r x and R and r 2 , <f> is the 
torsion angle. The diatomic reference Hamiltonians h ri (r x ) and h T2 (r 2 ) are defined as 


, , . h 2 d 2 ( . 
hri(ri) ~ + v(ri) ’ 

hr 2 ( r 2) = ~ — + V ( T ?)- 


(2) 


2^2 dr% 

Where V(r x ) and V(r 2 ) are the one-dimensional reference potentials for H 2 and CX and the 
Hi and ji 2 are the corresponding reduced masses. The ID potentials are the potential cuts in 
the respective coordinates with the system in the reactant chanel and all other coordinates 
are fixed at their equilibium geometry. 

The time-dependent wave function can be expanded in terms of the body- fixed (BF) 
rovibrational eigenfunctions defined in terms of the reactant Jacobi coordinates as 


E r2)Y/V‘(e I ,6 2 ,t), (3) 

TlUi Vo 3132 K 

where n in the translational basis label, and sine DVR basis are used as the translational 
basis, e is the parity of the system. <f> Ul and <f>„ 2 are the eigenfunctions of the reference 
Hamiltonia in Eq.(2). The angular momentum basis function Y ] J ^ i( {9 x ,9 2 ,(p) is the BF, 
parity adapted, total angular momentum eigenfunction, which can be written as 

Y X‘ = /,/ * x M nuK + £ (-l (4) 

3/2(1 -I- dtfoj 

where 


y,{iln,K = D J UK (a,M)y^ K (5) 

D J MK (a, (3, 7) is the normalized Wigner rotation maxtrix and is the angular momentum 
eigenfunction of j X2 
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with yjm the spherical harmonics. 

The split-operator method is employed to propogate the wave packet 

^ JMe (R,ri,r 2 ,< + A) = e~ iHoA/2 e~ iUA e~ iH ° A /2 'tp JMt (R ) r 1 ,T 2 ,t). 

The reference Hamiltonian H 0 is defined as 

h 2 d 2 


(7) 


H o 


0 2 ndR 2 

and the effictive potential operator in Eq. (7) 


+ h-n ( r i) + h r2 (r 2 ), 


( 8 ) 
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The initial wave-packet is chosen as the product of specific initial rovibrational states of 
the diatom H 2 and pesudo diatom CX and a localized standard Gaussian function <&k 0 (R) 


'k nol'ojo-ft'o — 9) — (77) dy 10 j 10 ( r 1 )d[/ 2n j 20 ( r 2 ) 1>: 'jlOj 2 ojl 2 oKo 


( 10 ) 


Fianlly, the initial-state-selected reaction probability is extracted from the time- 
dependent reactive flux in the standard way. 


B. Cumulative reaction probability and thermal rate constant 

In order to obtain an approximate full dimensional(12D) cumulative reaction probabil- 
ity (CRP) for J=0 and thermal rate constant, we first calculate the reduced dimensional (6D) 
CRP(J=0), denoted N J=Q ’ rd (E), by summing over all the initial state-selected rovibrational 
reaction probability P^ jlj2 C®) ’ 

N^(E)= £ < n > 

,^2=0 jij 2 =0 

where E is the total Energy. 


4 



Our reduced dynamics calculation has six degrees of freedom, while the H2+CH3 reaction 
system has twelve degrees of freedom. So the energy shifting method has to be applied to 
obtain the CRP for the full dimensionality using our reduced dimensional calculation. Thus 

N J =°’>““= N J =°' a (E - „). (12) 

here iqi/ 2 , ..., vq are the quantum numbers of the six vibrational degrees of freedom of the 
transition state, which are not included in our reduced dimensional calculation. These 
energies are given in the harmonic approximation by 

& + 0 3 ) 

i= 1 

The J-K shifting method is employed here to calculate the CRP for nonzero total angular 
momentum J. The total CRP, N(E) is defined as the sum of CPRs of all the open J and K 
channels 

j 

N(E) = + 1) Y, A - E'jk), (14) 

J - 0 K=-J 

where Ej K is the rigid rotor rotational energy of the reaction system at the transition state. 
In the present case, the system at transition state is a well defined symmetric top 

E x jk = B X J{J + 1 ) + (A t - B X )K 2 . (15) 


and B * are the rotational constants of H 2 CH 3 at the full dimensional transition state. 
The thermal rate constant can be computed exactly from the expression 

k(T) = / N(E)e~ E/kBT dE. (16) 

ft'UJreact J 0 

The above equation simplifies in the J-K shifting approximation, Eq. (14), to 

fit poo 

k(T) = - / N(E)e~ E/kBT dE, (17) 

h^react J 0 

where Ql ot is the rotational partition function of the reaction system at the transition state 
and Qreact is the reactant partition function, which is written as a product of vibrational, 



rotational, and translational partition functions. As noted in Ref. [24], the consistent ap- 
proach to calculate the vibrational partition function of CH 2 is to combine the results of the 
reduced dimensional quantum calculations with the ’missing’ vibrational frequencies which 
are not included in the reduced dimensonality. The purpose of doing this is to be consis- 
tent with the procedure that we use to convert the full dimensional CRP from the reduced 
dimensional CRP in Eq. (12). 

III. RESULTS AND DISCUSSION 
A. Numerical aspects 

The numerical parameters used in this work are: 60 sine basis functions for the transla- 
tional grid of R in the range of [2.7, 10.0] Bohr, 30 H 2 optimized vibrational basis functions 
for the T\ coordinates, 4 optimized CX vibrational basis functions ofr the r 2 coordinates, and 
total 1092 coupled angular momentum basis used for the angular part. The wave-packet was 
propogated with a time step of 15 for a total time of 6000 a.u. The rotational constants, A * 
and B t in Eq. (15), are 5.12 cm -1 and 2.00 cm -1 . The vibrational frequencies [24] that are 
’missing’ in our reduced dimensional calculations at the transition state are: 589 cm 1 (2), 
1439 cm -1 (2), 2959cm -1 and 3092cm -1 (2). These ’missing’ vibrational frequencies are used 
to obtain the full dimensionality CRP for J=0 in Eq,(12). 

B. Initial-state-selected probability 

Fig. 1 shows the reaction probability for the first four initial vibrational states of H 2 (vu 
ji=0) with CX at ground sta.te(v 2 =j 2 =0) as a function of kinetic energy. As seen, vibrational 
excitation of H 2 enhances this reactivity by lowing the reaction threshold. Fig. 1 also shows 
resonance-like structures in these vibrational excitation probability. And this feature has 
also been seen in the forward reaction, H + CH A [21,24], which indicates that the dynamical 
resonances exist in these reactions. The reaction probabitily of vibrational exciation of 
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CX(u 2 , j 2 = 0, i/ 1 =j 2 =0) are plotted in Fig. 2. Comparing to Fig.l, The reaction probability 
decreases as the vibrational excitation of CX increases. This is understandable due to that 
CX mainly funtions as a spector in this reaction. The more of the CH 3 vibrates, the harder 
for H atom in H 2 to ’catch’ the C atom in CH 3 . 

The rotational excitation probability of H 2 , with U\—l, u 2 =j 2 = 0 , for j\ up to 5 is given in 
Fig. 3 as a function of kinetic energy. In Fig. 3, ji=0 rotational state of H 2 gives the largest 
reaction probability. Comparing to Fig. 4, which is the reaction probability for different 
CX rotational states , j 2 =l initial rotational state has the largest reaction probability. This 
is because H 2 is much more vibrationally active than CX, and as a result, the isotropic 
rotational function of H 2 { iq=l, ji=0) has a larger radius, which ensures a bigger overlap 
between the rotational wavefuction of H 2 (ji—0) and C X (i/ 2 =j 2 = 0) than that of H 2 (j i=l), 
whose rotational wavefunction only distributes along the rq coordinate. So the reaction is 
more favorable to H 2 (ji=0) initial state in Fig. 3. On the contrary, the initial CX(j 2 = 0) 
state has a smaller radius of the isotropic rotational wavefuction, while the CX(v 2 —\, j 2 =\) 
initial state, which distributes along the r 2 bond, has the bigger overlap with H 2 (vi=ji= 0 ) 
at the linear transition state with C atom facing either H atom in H 2 . This explains why 
the CX{j 2 =\ ) initial state has the larger reaction probability in Fig. 4. Also one can see 
that the rotational excitations of H 2 have the bigger effect on the reaction probability than 
those of CH 3 . 

C. Cumulative reaction probability and thermal rate constant 

A normal mode analysis is performed to determine the ’missing 5 modes in our reduced 
dimensional quantum calculation, as also seen in Ref. [24]. Then Eq.(12) is used to obtain 
the full dimensional CRP for J=0. For total scattering energy up to 2.1 eV, all the missing 
modes given in Section A. contribute to convert the full dimensional CRP from the six 
degrees of freedom quantum calculation. And the full dimensional converted CPR is showm 
in Fig. 5 as a function of the total energy, which includes the difference between our six 
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degrees of freedom ZPE and the full dimensional normal mode ZPE, 0.936 eV. Using the 
harmonic approximation of the ’missing’ degrees of freedom in C H 3 is consistent with using 
this approximation in the energy shifting expression in Eq. (12). 

The resulting rate constant, obtained using J — K shifting is showm in Fig. 6 along with 
several different experiments. This figure shows that the present six dof calculation is larger 
than the experiments, espcially in the low temperature range. This feature, theoretical 
rate constant larger than the experiment, is also shown in the foward reaction H + CH 4 
[19,20,22-25]. Thus present study of the reverse reaction of H 2 + CH 3 supports the earlier 
studies on H + CH 4 and that barrier on the Jordan-Gilbert potential is too low. 

IV. SUMMARY AND CONCLUSIONS 

A six degree-of-freedom, time-dependent wavepacket propogation method, has been per- 
formed for the reaction, H 2 + CH 3 -4 H + CH i} by treating the reaction system as a 
pesudo-diatom-diatom reaction system. The energy-shifting method has been employed to 
obtain the full dimensional cumulative reaction probability for J=0. And J — K shifting 
approach is applied to this CRP to obtain the thermal rate constant which is compared with 
different experiments. A potential surface due to Jordan and Gilbert, which has been used 
in previous quantum calculation of the forward reaction H + CH 4 -4 H 2 + CH 3 , was used 
in these calculations. 

Initial state selected reaction probability and the cumulative reaction probability were 
calculated for zero total angular momentum. The vibrational excitations of H 2 bond substan- 
tially enhance the reaction probability. The vibrational excitations of CH 3 cause difficulty 
for H in H 2 to ’capture’ C atom, thus it has smaller effect on the reaction probability. The 
j 2 =■ 1 initial states of CH 3 gives the largest reaction probability comparing to ji = 0 of H 2 , 
this feature needs to be further investigated in furture work. 

Energy shifting method was used to convert the reduced dimensional CRP for zero total 
angular momentum to a full dimensional one. And our rate constant is larger than the 
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experiment, an indication that the barrier on the Jordan-Gilbert potential is most likely too 
low. 
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Fig. 1 The initial state selected reaction probability for the reaction H 2 {v\, ji — 0) + 
CHZ{y 2 = j 2 = 0) — > H + CH A as a function of kinetic energy (J=0). 

Fig. 2 Same as Fig. 1 but for the vibrational excitations of CH${v 2 ,j2 = 0) with H 2 (u i = 
3 \ = 0 )- 

Fig. 3 The initial state selected reaction probability for the reaction H 2 (v\ = 0, j\) + 
CH3(v 2 ,j 2 = 0) — >• H + CHi as a function of kinetic energy (J=0). 

Fig. 4 Same as Fig. 3 but for the rotational excitations of CHz{v 2 = 0,j2) 

Fig. 5 The cumulative reaction probability N J=0 ^ ull (E) as a function of total energy. 

Fig. 6 Comparison of the present reduced dimensionality thermal rate constant with different 
experiments. 1 Ref. [12] , 2 Ref. [9] , 3 Ref. [7] , 4 Ref. [11]. 
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